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ABSTRACT 

Context. With hundreds of exoplanets detected, it is necessary to revisit giant planets accretion models to explain 
their mass distribution. In particular, formation of sub-jovian planets remains unclear, given the short timescale for 
the runaway accretion of massive atmospheres. However, gas needs to pass through a circum-planetary disc. If the 
latter has a low viscosity (as expected if planets form in "dead zones"), it might act as a bottleneck for gas accretion. 
Aims. We investigate what the minimum accretion rate is for a planet under the limit assumption that the circum- 
planetary disc is totally inviscid, and the transport of angular momentum occurs solely because of the gravitational 
perturbations from the star. 

Methods. To estimate the accretion rate, we present a steady-state model of an inviscid circum-planetary disc, with 
vertical gas inflow and external torque from the star. Hydrodynamical simulations of a circum-planetary disc were 
conducted in 2D, in a planetocentric frame, with the star as an external perturber in order to measure the torque 
exerted by the star on the disc. 

Results. The disc shows a two-armed spiral wave caused by stellar tides, propagating all the way in from the outer 
edge of the disc towards the planet. The stellar torque is small and corresponds to a doubling time for a Jupiter 
mass planet of the order of 5 Myrs. Given the limit assumptions, this is clearly a lower bound of the real accretion 
rate. 

Conclusions. This result shows that gas accretion onto a giant planet can be regulated by circum-planetary discs. 
This suggests that the diversity of masses of extra-solar planets may be the result of different viscosities in these 
discs. 
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1. Introduction 

The detection of exoplanets is probably one of the most 
striking discoveries of the past 50 years in astrophysics. 
Up to now, more than 700 exoplanets have been found. 
This constantly increasing number allows a statistical 
approach to the analysis of exoplanets. The resulting 
statistical distributions are then crucial benchmarks for 
planetary formation models. For an exoplanet, the two 
easiest parameters to determine are the mass and the 
semi-major axis of its orbit. Low mass planets are still 
hard to detect, but it seems that the mass distribution 
of planets more than 10 Earth masses is double-peaked, 
with a peak close to the mass of Neptune, and another 
peak close to the mass of Jupiter Q. The first peak is 
probably an effect of observational biases as low mass 
planets are still hard to detect, b ut the second peak 
seems to be real (jMavor et al.ll20lil) . 

The most popular model fo r giant planet formation 
i s the core accretion model of Bodcnhe imer fc Pollack! 
(|1986l ). In this model, a solid core grows, embedded 
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in a gaseous proto-planetary disc . The gas within the 
Bondi radius (Bodcnhe imer et al.ll200fi l is gravitation- 
ally bound to the core, and forms an envelope. When the 
planet reaches a critical core mass of about 12 — 16 M® 
((Pollack et aBl996t iTarima fc Nakagawalll997D . the en- 
velope starts to contract quasi-statically and gas accre- 
tion rate increases, for an evolution timescale of a few 
million years. This long phase is followed by a runaway 
gas accretion stage occurring as soon as both core and 
envelope masses are approximately equal, during which 
quasi-static contraction and thus gas accretion rate dra- 
matically increase. This leads to an exponential growth 
of the planet, in a few 10 5 yr, which proceeds as long as 
there is gas available. 

This raises the question : what sets the terminal mass 
of a giant planet ? When a planet becomes massive 
enough (typically a Saturn ma ss), it starts to open a 
gap in the disc around i ts orbit dPapaloizou fc Linlll9"84t 
iLin fc Papaloizoull986t ICrida et al.ll2 006) . However, gap 
opening does not neces sarily imply that the gasflow onto 
the planet has ceased (jArtymowicz fc Lubowlll99"6l) be- 
cause the gap is not totally empty. Hydro-dynamical 



G. Rivier et at: Circum-Planetary discs as bottlenecks for gas accretion onto giant planets 



simul ations (iBrvden et al. | ll999t|Klevl l999t iLubow et all 
119991 : ILubow fe D'AngeldWoeJ) show that a Jupiter- 
mass planet keeps accreting gas almost as fast as if 
gap-opening had not occurred. Accretion stalls by gap- 
opening only once a mass in the range of 5 — 10 Jupiter 
masses is achieved. This is at odds with the masses of 
the giant planets in the solar system and of many of the 
extra-solar giant planets. In particular, the existence of 
Saturn-mass planets is a mystery in this scenario. 

An often proposed solution (e.g. : iThommes et al.l 

2008) is that the giant planets form in "dead zones" - 
regions in the disc with very low viscosity - and, con- 
sequently, open gaps that are much wider and cleaner 
than previously thought. These gaps could inhibit ac- 
cretion for planet masses of order of a Jupiter mass or 
less. However, this idea fails because, as demonstrated by 
ICrida et"aT1 ()2006[ ). the depth and width of a gap docs 
not depend just on viscosity but also on the temperature 
(i.e. scale height) of the disc. Thus, even in a dead zone 
the gap opened by a Jupiter-mass planet should not in- 
hibit significantly the accretion of gas into the planet's 
Hill sphere. Thus, we believe that the question on the 
terminal mass of giant planets is still open. 

We notice, however, that the iPollack et al.l (|1996( ) 
model considered an omnidirectional gas accretion to- 
wards the planet ; yet, it is physically impossible to reach 
this state due to angular momentum conservation of the 
inflowing gas. In fact, gas must form a consistent circum- 
planetary disc (CPD) once the mass of the planet is 
higher than about a hundred Earth masses and the plan- 
etary a tmosphere has shrunk well inside the planet's Hill 
radius (jWard fc Canudl2010D . This has been confirm ed 
by numerical simulations (e.g. lAvliffe fc Bate]|2009allbT ). 
Once a CPD is formed, the accretion rate of the planet 
depends on the ability of gas to lose angular momentum, 
either by its re-distribution within the disc (through vis- 
cosity) or exchange with external perturbers (e.g. the 
star). 

Most of the previous simulations considered a signif- 
icant viscosity in the CPD, comparable to that of the 
active zones of the circum-stellar disc. In this case, an- 
gular momentum is redistributed very quickly and the 
accretion of gas from the CPD onto the planet is ex- 
tremely fast, so that t he role of the disc can genera lly be 
neglected as shown by iPapaloizou fe Nelson! (120051) who 
assumed a viscosity about a — 10 -3 . However, it is pos- 
sib le that the CPD is i n an MRI dead state as suggested 
bv lTurner et al.l (|2010f ). leading to vanishing turbulence 
and a very low viscosity. Indeed, if the circum-planetary 
disc has a very low viscosity, then the transport of angu- 
lar momentum through this disc can be very inefficient 
and gas can only accrete onto the planet at the rate 
allowed by the removal of angular momentum by exter- 
nal perturbations, such as the stellar tide. Consequently, 
the circum-planetary disc may act as a bottleneck for 
gas accretion and dramatically slow down the planet's 
growth. If, as a result, the accretion of a planet does not 
enter the runaway phase but proceeds at a more regu- 
lar pace on a timescale comparable t o the circum-stellar 
disc lifetime (a few mi llion years) (Ha isch et al.l l200ll 
iHillenbrand et al.l 12008! ) then the observed mass spec- 
trum of the giant planets may be the consequence of the 
competition between gas accretion and gas dissipation. 



To test whether the idea of a slow accretion rate 
through a low-viscosity disc is realistic, we study in this 
paper the effect of the stellar tides. We are aware that 
Reynolds stresses from waves driven from the circum- 
solar disc may allow angular momentum losses in the 
CPD as well. However, the accretion rate due solely to 
the solar tide provides a lower-bound to the real accre- 
tion rate for a CPD of vanishing viscosity. If the result is 
encouraging — i.e. the mass-doubling time for a Jupiter- 
mass planet from solar tides is longer than the proto- 
planetary disc lifetime — then the idea of a low- viscosity 
CPD as a regulator of the gas accretion rate onto a giant 
planet is promising. In this case, future work will have 
to evaluate in detail the gas accretion rate with realistic 
3D hydro-dynamical simulations. Thus, our present pa- 
per should be considered as a first step in a long research 
plan. 

The structure of the paper is the following. In Sect. [5] 
we elaborate an idealised semi-analytical model to eval- 
uate the steady state accretion rate of a planet that is 
surrounded by an inviscid disc which undergoes verti- 
cal gas inflow and is submitted to a tidal torque from a 
distant star. Then we evaluate the two key parameters 
that enter in this model. The vertical gas influx rate is 
estimated in Sect. [3] from 3D simulations available in the 
literature. To compute the stellar torque, in Sect. 2] we 
perform 2D hydrodynamical simulations of a disc cen- 
tred on a Jupiter-mass planet. The disc feels the gravi- 
tational perturbation from the star, assumed to be on a 
distant, circular orbit. We measure the torque due to the 
stellar tides, after a steady-state is reached. With these 
two values in hand, in Sect. [5] we estimate the planet ac- 
cretion rate, expressed as a timescale for the doubling of 
the mass of a Jupiter-planet. Finally, we interpret these 
results and discuss their possible implications for giant 
planet formation theories in Sect. [5] 

2. Model 

The tidal torque and the accretion rate are not related 
in a trivial way in an inviscid disc. In this section, we 
present a simple model to derive an accretion rate on 
the planet from the tidal torque that we will measure in 
Sect. E 

2.1. Depletion of the CPD 

The torque exerted by the star on a ring of the disc is 
not entirely deposited locally in the disc. Only a fraction 
of the torque is deposited. This fraction depends on the 
disc viscosity and tends to zero for vanishing viscosity 
([Martin fc Lubowl l20lH ). The rest (i.e. the totality of 
the torque in the case of an inviscid disc) is passed to 
the a djacent rings through pressure effects (jCrida et al.l 
120061 . Appendix C). If no torque is deposited on a ring, 
then the gas in the ring does not lose angular momentum 
and does not move towards the central planet. 

It would be incorrect, though, to expect that the 
wave raised by the star does not promote gas accretion 
from an inviscid disc. The torque is transmitted from 
each annulus of the disc to its internal neighbour by 
pressure all the way to the inner edge of the disc (at the 
boundary with the atmosphere). This inner annulus, not 
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being supported by anything interior to it on a Keplerian 
motion, then loses angular momentum and falls onto the 
planet. Once the innermost annulus is emptied, the same 
fate occurs to the second innermost annulus and so forth. 
Thus, an inviscid disc is depleted from the inside out as 
gas is accreted onto the planet's atmosphere. 

More precisely, let us note X the radius of the region 
of the CPD that is emptied by the torque over a time t. 
The total orbital angular momentum of the gas present 
within r < X is given by : 



mass flow inside r < Xc ■ 



L = f S27T7- x r(rfl) dr 
= —V^/GM~ P X 5 / 2 





(1) 

(2) 



(assuming the surface density of the gas, E is indepen- 
dent of r). Here, Q — J is the angular speed of the 
gas in Keplerian motion at a radius r, with M p the mass 
of the planet. 

By definition of X, L should be exactly opposite to 
the total angular momentum taken by the star, which is 
the torque T, accumulated over time : T x t. This gives : 



X(t) = 



\T\ 



2/5 



.2/5 



(3) 



We see that the radius X increases with time, as ex- 
pected. 

2.2. Infilling of the CPD 

Here, we assume that the CPD is non-viscous. Thus, 
the CPD doesn't spread radially, and a ny gas input at 
i ts ou ter e dge just stays th e re. Ho wever, iMachida et al.l 
(2008) and iTanigawa et alJ (|201l| ) have shown that ver- 
tical inflow onto the CPD is the main source mass of the 
CPD. Let us call C this vertical mass flux per surface 
unit of the CPD. 

Then, the surface density in an initially void region 
is £ = Ct. Input into Eq. ©, this gives: 



X c 



ITI 



2/5 



(4) 



This implies that there is a radius Xc (independent 
of time) such that inside r < Xc, depletion due to the 
stellar torque T and infilling due to the vertical inflow 
C exactly balance out. In an arbitrary time span, all 
the gas input within a circle of radius Xc around the 
planet has exactly the orbital angular momentum taken 
by the stellar torque. As a consequence, the whole gas of 
this region falls onto the planet while the vertical inflow 
completely refills it in the same time. Thus, the region 
within Xc reaches a steady-state density. 

The disc outside Xc, however, does not accrete and 
gas piles up. This issue will be discussed later in Sect. 12.41 

2.3. Accretion rate estimate 

We can now write the planet accretion rate as a function 
of the torque. The mass flow on the planet is simply the 



CnXr 



4/5 



M„ = 7T 



(5) 



Note that the accretion rate depends weakly on the gas 
injection flux C, whereas it depends almost linearly on 
the torque exerted by the star T. This equation is cen- 
tral in this paper, as it links the stellar torque to the 
planetary accretion rate. 

2.4. The issue of gas-piling beyond Xc 



The model introduced in Scct. l2.2l raises the issue of gas- 
piling in the outer part of the disc. Indeed, if the density 
in the outer disc continued to increase, the total torque 
T exerted through the wave would increase as well, as 
it is proportional to E. This would imply an increase of 
the accretion radius X , and of the accretion rate M p . 
Moreover, the density in the outer disc could reach the 
gravitational instability limit, changing completely the 
structure of the CPD and leading to outburst s of accre- 
tion, as suggested bv lMartin fc Lubowl (|201ll ). 

In reality, though, gas cannot pile-up indefinitely. If 
the density in the CPD becomes too large, gas is repelled 
by pressure and the inflow from the circum-stellar disc 
has to stop. Thus, in this section we estimate the density 
in the outer CPD, by imposing the pressure equilibrium 
between the CPD and the circum-stellar disc (CSD). Our 
approach is detailed below. 

We assume the CSD is much thicker than the CPD. 
This, and the fact that the inflow fr om the CSD to 
the C PD is along the vertical direction ([Tanigaw a et al.l 
1201 ll ) suggests that the CSD "weights" on top of the 
surface of the CPD. Then, at equilibrium, the pressure 
in the CPD at a height PTcpd has to be equal to the 
pressure of the CSD at the same height. 

Let us now assume a standard Gaussian pressure 
profile in both discs: P(z) — P exp (— \{jj) 2 ) where 
Pq is the pressure at the mid-plane (denoted hereafter 
P CSD and P CPD for the CSD and CPD respectively). 
At height H C pd, the pressure in the CSD is P CSD , 
given our assumption that Hqsd ^> Pcpd- In the CPD, 
P{H C pd) = P^ PD /Ve- p o and E are linked by the equa- 
tion of state Po = CgE (where c s = H£l is the sound 
speed). The equation Pcpd(Pcpd) = Pcsd(0) can then 
be solved for the sole unknown of the problem : 



E C pd M = ve -77771 



M Q r 



M p r Q 



— Ecsd 



where Mq is the mass of the star and r Q the semi-major 
axis of the planet. 

For (H/r P) ) n Sr > = 0.05 dPietu et alJ I2007D . 
(P/V)cpd = 0.3 (|Avliffe fc Batell2009bD M^~ = 10 3 M P , 
r = 14.4 Rh (see simulation set-up in I4.1T[) and r = 
0.3 Rh which corresponds to the truncation radius of the 
CPD, as explained in 14.1.31 one gets Pqsd = 8 x PTcpd, 
justifying our assumption. This gives : 



^cpd 



(r) = 



0.31 R H 



^CSD 



(6) 
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with Rh the Hill radius of the planet : Rh = 
r MM' 73 

So, we estimate that the surface density of the CPD 
at 0.3 Rh is about equal to the local surface density of 
the CSD, and can not exceed this value. This value of the 
density will be used in our simulations for tidal torque 
calculation in Sect 14. 1.21 

3. Gas inflow estimate 

To apply the method outlined in Sect. [2j we need to 
estimate C: the flux per unit area of the gas injected 
onto the CPD . To do that, we use the results from 
lAvliffe fc Bate] (|2009bD . Those simulations give the flux 
of gas through the Hill sphere (Mjj, identified in that pa- 
per as the "planet accretion rate" ) . For a planet of mass 
M p = 333 M© (which is consistent with o ur choice of pa- 
rameters) and a locally isothermal disc, Ayliffe & Bate 
(|2009bD find M H = 8 x 1C T 5 M, IuD /yr. Given tha t almost 
all of this flux is vertical (jTanigawa et al.ll201l[) . we can 
set Mr — CnR H , which gives 

C = 2.54 x 1(T 5 M Jup /yr/R H 2 . (7) 

4. Numerical simulations and Torque 
computation 

The model introduced in Sect, fallows us to derive an- 
alytically an accretion rate onto the planet, given a gas 
inflow on the whole disc and an accurate knowledge of 
tidal effects from the star, i.e. an estimate of the stellar 
torque on the disc. Having evaluated the former in the 
previous section, here we describe how we compute the 
torque numerically. We start by reviewing the set-up of 
our simulations. 

4.1. Simulation set-up 
4.1.1. Code 

To run simulations, we have used the code FARGO0 
from lMassetl (|2000alibh . which is a 2D Eulerian, non self- 
gravitating, isothermal code, with a polar grid. FARGO 
is very suitable for simulations of the interactions be- 
tween planets and the proto-planetary disc. Here, we 
focus on the environment of the planet. Thus, we have 
slightly modified FARGO to adapt it to our case of study. 

From now on, the central body will be the planet, 
while the star will be the orbiting body around the planet 
and its disc. Simulations take place in a planetocentric 
frame, corotating with the star. The orbit of the star is 
circular and fixed. 

The mass unit is the mass of the central body : the 
planet, M p . The distance unit is Rh- The time unit is 
where Hh — \JGM P /R H is the angular speed at 
r = Rh- Note that the star's angular speed is thus fi© ~ 

/ GMq _ Q H 
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The mass of the star is taken to be M© = 1000 M p , 
which makes our system scaled like the Jupiter-Sun sys- 
tem. This gives r Q = 14.4 Rh- 

The star being out of the grid, we do not apply any 
smoothing in the computation of the force exerted by 
the star onto the disc. This is actually necessary for an 
accurate computation of the differential, tidal forces. 

4.1.2. Gas parameters 

In our simulations, we take an initial density profile 
£(r) = So (^) with E = 5 x 10" 4 M P /R 2 H . 

This is in agre e ment w ith the global simulations from 
lAvliffe fe~ Bate (2009a). Moreover, the surface density 
evaluated at 0.3R H is ~ 9x 10~ 4 M p /R 2 H as a — 1/2 (see 
below). This should be compared to the surface density 
at the location of Jupiter: Ecsd ~ 10~ 3 Mj uu /R% in 
the Minimum Mass Solar Nebula (|Weidcnschillingj |l977t 
lHavashi|[l98l . Thus, our disc fulfils the condition of 
maximum pile-up, that we evaluated in Sect. 12.41 

Due to numerical reasons, FARGO does not allow 
the modelling of a perfectly inviscid disc. Thus, we take 
a low but non-zero kinematic viscosity v = 10~ 5 R H flH, 
unless specified otherwise and adopt a surface density 
profile a — 1/2 so that the viscous torque between ad- 
jacent annuli is zero. Notice, moreover, that our method 
to compute the gas accretion onto the planet does not 
use the actual flow of the gas, but only the shape of 
the spiral wave generated by the star (whose structure 
is detailed in Sect. 15. 1[) . Thus, our results should not be 
significantly affected by our choice of as we will check 
in Sect. O 

The Equation Of State in FARGO is locally isother- 
mal, with P = c s 2 E, where c s is the sound speed and 
c s = HVt. The aspect ratio, H/r, is constant in all of our 
simulations and independent of r. 

Th e H/r of the CPD has been measured to be 0.3 — 
0.4 in lAvliffe fc Bate! (|2009al) . Unfortunately, the code 
FARGO is not stable for values of the scale height that 
are so large. Thus, we have run simulations with H/r 
values ranging from 0.05 to 0.2, in order to extrapolate 
the results to H/r = 0.4. Fortunately, as we will see 
below, the cumulated torque at the inner edge of the 
disc turns out to be independent of the H/r value, which 
makes the extrapolation trivial. 

4.1.3. Numerical parameters 

In several previous studies, simulations have shown that 
the disc reaches an equilibrium distribution with a sha rp 
truncation at its oute r edge (jAvliffe fc Bate! l2009al lbl; 
iMartin fc Lubowll2011l) . This truncation radius is about 
0.3 — 0.4 Rh a nd is likely set by tidal t runcation effects 
as outlined by IMartin fc Lubowl (|201ll) . Indeed, we ob- 
served ourselves such a truncation radius in simulations 
with a grid extended to the Hill radius. 

Outside the truncation radius, the motion of the 
gas relativ e to the planet i s stron gly non-Keplerian, as 
noticed by lAvliffe fc Bate! (|2009al) . The interaction be- 
tween this region and the CPD is not clearly under- 
stood. Therefore, as we just want to measure the effect 
of stellar tides on the CPD, we consider in this paper 
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a limit case where the CPD is completely isolated and 
we truncate the grid at i? ut = 0.3 Rh. The inner bor- 
der radius is R in = 0.01 Rh ; for a Jupiter mass planet 
at 5.2 A.U., it corresponds to about 540 000 kilometres, 
which is smaller than the semi-major axis of Europa. 

We choose a logarithmic radial spacing for the grid 
(the width of a ring 5r is proportional to its radius r), 
which is suitable to have an accurate resolution very 
close to the planet. The resolution is taken such that 
59 = Sr/r < Hjhr (see Tabled]). This resolution is ap- 
propriate to resolve the pressure density wave exerted by 
the star, as H is the typical pressure scale length. 

4.1.4. Boundary conditions 

We use a non-reflecting boundary condition for both the 
inner and the outer borders of the disc. This prevents 
density waves from reflecting at the border of the grid. 
It is important for an accurate evaluation of the torque 
onto the disc through the wave. It also does not allow gas 
inflow or outflow and keeps the mean density on the first 
and last rings constant with time, which is consistent 
with our model of an isolated disc. 

4.2. Torque computation 

The gravitational torque exerted by th e star on each ring 
of the CPD is computed as follows (see lMartin fc Lubowl 
|2011D . First, it is computed on each grid-cell. To the 
direct torque T g , we have to subtract an indirect torque 
Ti due to the acceleration of the frame (i.e. of the planet). 

Denoting by r — (x c , y c ) the coordinates of a cell of 
mass m ce \i, by r Q = (xq, y@) the coordinates of the star 
and by d = r — r Q = [xd, yd) the mutual distance vector, 
the direct and indirect forces exerted by the star on the 
cell write: 

~ ^ M\ 

_ GM m ccn r Q 



Thus, the gravitational torque exerted from the star on 
one cell is : 

imi = \\T g -Ti\\ 

= ||rA(/ ff -/i)|| 

= -[XcUg-Vd ~ fi-Vo) ~ Vc(fg-Xd - fi-X®)] 

Then, we compute the azimuthal sum over all the 
cells of an annulus, to obtain the gravitational torque 
felt by one ring. Finally, the total torque on the disc 
is the sum of the torques felt by the individual rings. 
Thus, it is useful to introduce the notion of cumulative 
torque, which is the sum of the torques felt from the 
first ring to the ring in consideration. The total torque 
is therefore the cumulative torque at the last ring. For 
our purposes in this paper we consider the outermost 
ring to be the first and the innermost one the last (un- 
like iMartin fe Lubowl 1201 lL who adopted the opposite 
convention) . 




-0.2 -0.1 0.1 0.2 



Fig. 1: Density map for an aspect ratio of 0.15 ; the color 
corresponds to the logarithm of T,/(Mj up R H 2 ). The disc 
is extended from r = 0.01 to 0.3 Rh- 



5. Results 

In this section, we describe the results of numerical sim- 
ulations on the structure of the CPD and present our 
estimate of the stellar torque. Then, we derive results 
for the accretion rate according to our analytical model. 

5.1. Disc structure 

As expected, we see a two-armed spiral density wave, 
created by the tides of the star. The wave propagates 
radially inwards with the speed of sound. As c s = Hi}, 
its radial velocity is H/r times its azimuthal velocity. 
Thus, the wave has the shape of a logarithmic spiral, 
with its pitch angle given by H/r. This is clearly seen in 

Fig.m 

5.2. Torque 

We measure the torques after the simulations have 
reached a steady - state (after 150 orbits at 0.3 Rh)- 
IMartin fc Lubowl (|2011[ ) showed that the total torque 
is negative and that most of it is exerted on the outer 
region of the disc (which represents less than 10 % of 
the disc mass) . Therefore, the variations of the cumula- 
tive torque in the inner part of the disc (r < 0.1 Rh) 
are negligible , as can be seen on Fig. [2j This figure 
shows the cumulative torques for various H/ r, after 250 
orbits at r = 0.3 Rh- T hey oscillate from r = .3 to 
0.1 R H - As explained by IMartin fc Lubowl (|201l[ ), the 
disc can be divided in four quadrants, starting from the 
line connecting the planet to the star, rotating in the 
anti-clockwise direction. When the two-arm wave is in 
the first and third radiant, the torque is negative; when 
it is in the second and fourth one, the torque is positive. 
Thus, the oscillations of the torque correspond to the 
passage of the wave from even to odd quadrants and vice 
versa, due to its spiral shape. The wave-length of these 
pseudo-periodic oscillations is the radial propagation of 
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Fig. 2: Cumulated torques (in Mj up R H il H ) as a function 
of radius (in Rr) for different aspect ratios. 



the wave during half a revolution around the planet, that 
is 7r r (H/r). 

Notice that, although the radial profile of the cumu- 
lative torque is different from simulation to simulation 
given that the wavelength of its oscillation is a function 
of H/r, the limit value at the inner edge of the disc is re- 
markably insensitive on H/r. Thus we can safely assume 
that the total cumulative torque for a H/r = 0.3 — 0.4 
disc (the realistic valu e of the scale height according to 
lAvliffe fc Bate1[2009al is also - 4 x 10" 8 . 

We can now check whether our choice of the disc vis- 
cosity v impacts significantly the total torque that we 
measure. For this purpose, Table [1] shows the torques 
measured in simulations with various viscosities. Some 
differences are visible, which is a sign that our simula- 
tions are not dominated by numerical viscosity. However, 
the total torques are very similar to each other (within 
10% for an order of magnitude change in v). This gives 
confidence that the total torque that we estimate is valid 
also in the limit of an inviscid disc. 



5.3. Accretion rate 

In Sect. [31 we obtained a value for the influx of gas onto 
the CPD C. Then, thanks to numerical simulations, we 
measured the stellar torque for different aspect ratios 
and low viscosities and obtained its order of magnitude 
in Sect. 15.21 Given these values, we can derive analyt- 
ically the planet accretion rate M from Eq. ([5]). The 
resulting values are listed in Table [1] along with results 
on the torques and the corresponding steady-state ra- 
dius Xc- We find that Xc is of the order of ~ 0.05 Rr , 
approximately 50 times the current physical radius of 
Jupiter. 

As can be seen in Table [TJ the accretion rate hardly 
depends on the aspect ratio, and only increases by about 
5% when H/r doubles. This is because, as we said above, 
the total cumulative torque is insensitive to the H/r 
value of the disc. The viscosity of the gas in the sim- 
ulation also has very little influence on the measured 
torque. In all the cases, we find that the accretion rate 



This accretion rate is low, making Jupiter double its 
mass in 5 million years. The torque that we measure is 
of course proportional to the surface density of the CPD 
in the simulations, but we have shown that the adopted 
value is realistic, and could hardly be larger. The param- 
eter C may be poorly constrained, but the dependence of 
M p on C is such that two orders of magnitude difference 
in C would only change M p by a factor 3. Therefore, 
this order-of-magnitude estimate of the accretion rate in 
an inviscid CPD solely perturbed by the star, is robust. 
As a consequence, the accretion of giant planets could be 
much slower than expected, thus preventing the runaway 
accretion phase of the planet. 

6. Conclusion 

The classic mode l for giant planet formation 
(|Pollack et al.l Il996| ) predicts that the final phase 
of gas accretion occurs in a very fast runaway mode. 
Consequently, planets should keep accreting gas until 
they are so massive they open very wide gaps in the disc, 
which occurs when they reac h a mass equa l to m u ltiple 
times the mass of Jup i ter (iBrvden et all 119991; iKlevI 
[T^lLubow et al.lll999HLubow fc D'Angeloll200CT . 

In this paper we considered that, during the alleged 
runaway growth phase, a planet should be surrounded 
by a circum-planetary disc (CPD) due to angular mo- 
mentum conservation of the gas global flow. Thus, most 
of the gas that it accretes should have passed through 
this disc. The viscosity in the CPD may be very low, 
if the planet is located in a dead zone. Therefore, we 
have investigated whether a non-viscous CPD could act 
as a bottleneck for gas accretion. If, consequently, the 
gas accretion timescale can become comparable to the 
timescale of gas removal from the proto-planetary neb- 
ula, the observed large spread in giant planet masses 
could stem from the competitions between these two 
timescales. 

Considering that there is no radial drift in a non- 
viscous disc, we have developed a model for a steady- 
state non- viscous CPD. The disc is fed by a vertical gas 
inflow from the surrounding environment (as previously 
observed in 3D numerical simulations, see Machid a et al.l 
[20081: iTanieawaet "all [2011 lAvliffe fc Bate! I2012D . un- 
til a pressure equilibrium is reached. The surface den- 
sity of the CPD is thus analytically determined and 
found in agreement wi th previous numerical simulations 
(lAvliffe fc Batdl2009ah . 

However, the CPD is perturbed by the star. This re- 
sults in the formation of a two-armed logarithmic spiral 
density wave, that propagates all the way inwards down 
to the inner edge of our disc. As a consequence, a nega- 
tive torque from the star is deposited in the very inner 
regions of the disc, where gas consequently falls onto the 
planet. We find that the planetary accretion rate de- 
pends almost linearly on the cumulated stellar torque, 
and weakly on the vertical gas inflow (Eq. ([S]) ). 

Running 2D simulations of an isolated disc, in a plan- 
etocentric frame, extended to 0.3 Hill radius, we have 
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Table 1: Simulation parameters and results. N r and N s are the number of radial and azimuthal grid-cells respectively. 
v is the viscosity, T the cumulated stellar torque, Xc the equilibrium radius between the inner part of the disc in 
steady-state and the outer part. Finally, M p is the accretion rate onto the planet. 



studied the effect of the star on the disc, and measured 
the torque. We find that the torque is negative and basi- 
cally independent on the disc's aspect ratio. This allowed 
us to derive the accretion rate of a giant planet sur- 
rounded by a non-viscous CPD perturbed by the star. 
We find that the mass doubling time for a Jupiter-mass 
planet is about 5 Myr. This timescale is much longer 
than t hat in the runaway accretion mode of lPollack et al.l 
(1996). However, such a low accretion rate is valid in 
the limit condition of an inviscid disc. In reality, even 
in a dead zone the viscosity is not null, though it is 
still extremely hard to estimate quantitatively its value. 
Furthermore, viscosity, if related to ionisation, can in- 
crease with time as the disc becomes less optically thick 
(|Turner et alj|2010ft . As a consequence, the broad mass 
range of observed exoplanets may come from a range 
of viscosities (or viscosity evolutions) in their original 
CPDs. 

Moreover, we have only studied the effect of the star 
on gas accretion inside the CPD. Several other effects 
may intervene to make the CPD lose angular momen- 
tum, which need to be investigated in the near future. In 
particular, interaction with the gas beyond 0.3 Rh from 
the planet and outside the Hill sphere may perturb the 
flow in the CPD. The study of such interactions would 
require further investigations in global simulations and 
is beyond the scope of this paper. 

In conclusion, it emerges from this paper that the 
accretion history for planets more massi ve than Saturn 
(the m ass beyond which a CPD forms |Ayliffc fc Bate! 
2009b) may be dominated by the viscous evolution of the 
CPD. We speculate that planets with masses between 
Saturn's and a few times Jupiter's may have formed 
in "dead zones" with different levels of low viscosity in 
their CPDs. Instead, very massive planets (5-10 Jupiter 
masses), which reached the mass limit of the runaway 
growth process imposed by gap opening in the disc, 
should have formed in active zones, where CPDs did not 
act as bottlenecks to accretion. 
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